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Abstract 

We describe an algorithm for the continuous monitoring of the biomass 
and ethanol concentrations as well as the growth rate in the Mezcal 
fermentation process. The algorithm performs its task having available 
only the on-line measurements of the redox potential. The procedure 
combines an artificial neural network (ANN) that relates the redox 
potential to the ethanol and biomass concentrations with a nonlinear- 
observer-based algorithm that uses the ANN biomass estimations to 
infer the growth rate of this fermentation process. The results show 
that the redox potential is a valuable indicator of the metabolic activity 
of the microorganisms during the Mezcal fermentation. In addition, 
the estimated growth rate can be considered as a direct evidence of 
the presence of mixed culture growth in the process. Usually, mixtures 
of microorganisms could be intuitively clear in this kind of processes, 
however the total biomass data do not provide definite evidence by 
themselves. In this paper, the detailed design of the software sensor 
is presented, as well as its experimental application at the laboratory 
level. 
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1 Introduction 



Mezcal is a Mexican distilled spirit produced from the core of the Agave 
plant (the "pihas"). Several species are used for the Mezcal production, i.e., 
Agave salmiana, potatorum, angustifolia, tequilana, etc, and each produces a 
slightly different Mezcal. Therefore, the tequila can be considered a regional 
type of Mezcal, restricted to the usage of Agave tequilana as raw material pQ. 
It is important to mention that during this Mezcal fermentative process the 
syrup obtained from the juice of cooked "pihas" is left to naturally (spon- 
taneously) ferment. By natural alcoholic fermentations we refer to the ones 
that start by themselves when a wide spectrum of different microorganisms 
starts fermenting. Therefore, the fermentation of Agave syrup into Mezcal 
is a complex biochemical process involving a whole realm of interactions 
between microorganisms. 

During fermentation, the microorganisms employ sugars and other con- 
stituents of Agave syrup as substrate for their growth, converting them into 
ethanol, carbon dioxide, higher alcohols and their esters, and other metabolic 
compounds that contribute to the chemical composition and sensory quali- 
ties of the Mezcal [2]. We also notice that the qualitative and quantitative 
composition of the microbiota in fermenting musts could depend on the fol- 
lowing factors: region of the Agave origin, production procedure, initial cell 
concentration, temperature, and ethanol concentration. As consequence, 
the organoleptic properties are also the result of the diversity and compo- 
sition of microorganisms and their dynamics and frequency of occurrence. 
Consequently, this fermentative process is a vital stage in Mezcal making. 
Thus, it is very important to know more about the dynamics of the en- 
tire microflora during the alcoholic fermentation process. In other words, 
an algorithm that enables a monitoring process could be fundamental for a 
quality control that ensures, at least, homogeneity in the final product. 

Nevertheless, a bottleneck in all biochemical monitoring process is often 
the lack of sensors for biological variables. Moreover, it is a well-known 
issue that in order to monitor many biotechnological processes, the prob- 
lem of growth rates estimations represent a strategic feature. That is why 
several techniques have been developed to estimate on line the biological 
variables from the available measurements, which are usually dealing with 
physico chemical variables. Depending on the obtainable information about 
the process, there exist many possible types of estimators that can be used 

0, m, 0. 

Previously, various attempts of relating the redox potential to fermen- 
tation processes have been made taking into account that redox potential 
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assesses the growth ability of microorganisms, as well as the physiological 
activity in a given environments [6], [7J, [H]> IS]- Particularly, the practical 
significance of redox potential and oxygen content at various stages of wine- 
making was examined by [TO]. Many chemical, enzymatic and biological 
processes in wine are correlated with the oxidative state of the wine. 

Monod kinetics, which were originally derived from laboratory experi- 
ments with pure cultures and single substrates, are frequently applied to 
describe the behavior of undefined mixed cultures growing with single sub- 
strates or complex substrate mixtures [TT], [TO], [TO]. In this case, the growth 
parameters that have been used represent overall values reflecting the growth 
constants of many different strains with respect to the multicomponent sub- 
strate and the frequencies and concentrations of both the different substrates 
and microbial strains. Nevertheless, it is well known from control processes 
that an accurate model leads to better control design and therefore to bet- 
ter closed-loop performance. There is recent progress shedding light on the 
dynamical processes underlying the growth of mixed culture in a mixture of 
substrates [TO], [TO]. However, under such complex situations, the accurate 
determination of the growth rates is still an open subject. 

The work is organized as follows. Section[2]is devoted to a concise presen- 
tation of the fermentation experiment performed to illustrate our approach. 
The software sensor, which is a combination of an ANN and an adaptive 
observer scheme, is described in subsection 12.41 The results obtained by 
applying this software sensor to the Mezcal fermentation process together 
with a number of comments are included in Section [3l The paper ends up 
with the concluding remarks. 

2 Materials and Methods 

2.1 Microorganism and culture conditions 

In order to evaluate experimentally the performance of the estimation algo- 
rithm, we performed six individual batch experiments using inocula of native 
microorganisms (without the addition of any commercial strain). The must 
(Agave syrup) was obtained from Agave salmiana, a species from the Mex- 
ican plateau (or altiplano) of the geographical region of San Luis Potosi. 
This must was centrifuged at 8000 rpm x 10 min and stored in a frozen 
state at — 20°C prior to experiments. 

The batch fermentations were carried out in a bioreactor (Applikon, 
Schiedam, the Netherlands) of 1 liter. The bioreactor is equipped with pH 
and redox sterilizable electrodes (Pt-Ring, Applisens, Sensor Innovation, 
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Applikon). The electrodes are connected to a console for data acquisition 
(Bioexpert, Data Acquisition Control Program, Version l.lx, Applicon), a 
device which is connected to a computer where the data are stored and 
computed. 

The bioreactor was filled with 900 mL of must as a culture medium, 
100 mL of the inoculum in its exponential growing phase (biomass 0.1 g/L) 
and 0.1 % of ammonium sulfate at final concentration. The initial condi- 
tions of the fermentation were settled at a temperature of 32.5°C and initial 
sugar concentration of 70 g/L. The pH does not show a dynamic evolution, 
maintaining itself at a value of 4 during the whole process. For a schematic 
representation of the experimental setup and the data acquisition system 
see Fig. [TJ 

2.2 Analytical procedures and measurements 

The batch processes have been monitored for 14 hours, through sampling 
under sterile conditions. In order to quantify biomass and ethanol concen- 
trations, 5 mL samples of culture was removed every 30 min. The samples 
were cleared by centrifugation at 6000 rpm for 5 minutes at room tempera- 
ture. The next step was to collect the supernatant phase and store it frozen 
at — 20° C prior to be analyzed. The obtained pellet was resuspended in 
distilled water in order to proceed with the biomass analysis. 

2.2.1 Biomass and ethanol quantification 

The biomass measurements have been performed using (varian) UV spec- 
troscopy at 600 nm. The obtained values were interpolated with a standard 
curve of cell dry weight concentration. 

For the determination of the ethanol, 1 mL of each sample with a final 
dilution 1:10 and 1 mL of 1-butanol in vortex motion for 5 min and followed 
by centrifugation at 6000 rpm for 5 min. The organic phase was analyzed 
by gas chromatograpy (6890N Agilent technologies, Wilmington, DE) pro- 
vided with a capilar column HP-Innowax (30 m x 0.25 mm i.d., 0.25 m film 
thickness; Agilent technologies, Wilmington, DE) and an auto-sampler 7863 
(Agilent technologies, Wilmington, DE) with a split relation of 25 : 1. 

The chromatographic conditions were 35°C for 2 min, increased at the 
rate of W°C/min up to 80°C, and maintained at the latter temperature for 
15 min. The carrier gas was helium at a flow rate of 1.5 mL/min. The 
temperatures of the injector and flame ionization detector (FID) were set at 
220° C and 250° C, respectively. The ethanol concentration of the samples 
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was determined by means of a calibration curve of known standard solutions 
of ethanol. 

Data acquisition device 
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Personal Computer 
with the Software Sensor 

Figure 1: The schematic representation of the experimental setup in our 
laboratory. 



2.3 Redox potential 

The measurement of redox potential is relatively fast, accurate, and reliable 
and its values give an insight into the oxidation process as well as the in- 
herent ability of reduction in the process, which is well established in the 
case of wine [ID]. We mentioned that improvements in manufacturing pro- 
cesses and equipment have allowed beer, wine, as well as Mezcal producers 
to increase both the quantity and quality of these products. It is crucial 
for the commercial competition to maintain good quality control practices. 
Higher good-quality productions of the manufacturing plants require faster 
and more sophisticated analytical techniques to achieve them. 

The measured values of the redox potential can give information on redox 
reactions in wine, which have an important influence on its quality and 
stability [10J. During storing and aging of wine, oxidation and reduction 
processes affect the character and taste of wine to a considerable extent 
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|10| . [9]. We surmise that the same features could occur in the case of the 
Mezcal making. Checking this assumption is one of the main goals of this 
research. In the case of tequila, there exist previous efforts to understand 
the relation of the organoleptic properties with the process stages [T7] , [18] , 
|19| . [20] . However, both for Tequila and Mezcal there is no previous report 
of an online-measured variable that allows its fermentation monitoring in 
real time. 

The redox potential measurements were acquired periodically each 0.01 hr 
during 14 hrs, and the data were stored and computed on-line in a PC (see 
Fig.E}. 

2.4 The Software Sensor 

By the software sensor we mean the algorithm generated by the coupling 
of the ANN and the adaptive observer. In this section we develop step by 
step this type of algorithm. First, we present the mathematical model that 
stays as the background of our approach. Then we explain the relationship 
with the measured output, the redox potential, and the ANN scheme. We 
end the section with a discussion of the adaptive scheme used to infer the 
growth rate function. 

2.4.1 Mathematical model of the fermentation process 

Batch microbial growth in a completely stirred bioreactor is commonly de- 
scribed by the following differential equations [3], [E] 

-k lf i{t)X 1 
k 2 n(t)Xx 

where X\ represent the total biomass concentration ([g/l]), X2 refers to the 
total substrate concentration ([#//]) and X3 is the product concentration 
([g/l]) (ethanol in this case), k\ represents the biomass yield, k 2 the product 
yield and [x (t) is the growth rate. All the parameters that have been used 
represent overall values reflecting the many different strains with respect to 
the multicomponent substrate. It is important to note that the scheme here 
proposed do not need the values of k\ or k 2 and /i(t) is inferred by means of 
the software sensor. 



dXi 

~dT 

dX 2 

dt 
dX 3 

dt 
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2.4.2 Artificial neural network 



In order to relate the redox measurements to the ethanol and biomass con- 
centration an ANN procedure is applied. The methodology that we carried 
out includes a forward-propagation training algorithm for the ANN using 
some of our experimental data. In order to perform our task we construct 
a model of the following form: 

Xi = /(Xt) , (1) 

where X4 represents the redox potential measurement data (mV), X\ is 
the set of biomass concentration data (mg/L) and the function /(A4) is 
approximated by means of the ANN procedure. The ANN architecture is of 
the standard type [22] with a single ANN hidden-layer containing 10 units. 
The same scheme was used for the case of the ethanol but using A3 = /(A4) 
as the ANN model. 

Each unit of this network uses a sigmoid function as the activation func- 
tion. On the other hand, the output contains a linear activation function, in 
our case the identity. The feed-forward training algorithm considered here 
is the conjugate gradient method [23]. Three of the six individual batch 
experiments were used to provide data for the training process. The error 
given by the ANN after the training was as low as 0.0029. 

2.4.3 The adaptive observer 

The analysis of the adaptive observer scheme used here is based on the 
following realistic assumptions: 

(Al) The specific growth rate (j,{t) is positive and bounded, that is an un- 
known Umax exists Such that < fJ,(t) < fJ-max- 

(A2) There is no growth without substrate: X2 = => fi(t) = 0. 

(A3) The time derivative of [i(t) is bounded: | d 4^ | < Mi, where Mi E R+. 

2.4.4 On-line estimation of fj,(t) from measurements of X\ 

Since X\ is available through the neural algorithm mentioned previously, we 
can now rewrite the output of the system as: 

y s = x u (2) 
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where y s is the set of on-line measurements of the system which is available 
indirectly through ANN means. Then, according to the discussion in [3], 
the following algorithm can be used to estimate fJ,{t)'. 



d Xi 



y s jl (t) + K-i y s [y s - X x 
K-2Vs [Vs ~ Xi) , 



) 



d /2 



dt 



where X\ and jl represent the estimated value of X\ and respectively. 
The constants K,\ and /C2 must be chosen such that: 



The above condition ensures the asymptotic convergence of the observer 
error to a neighborhood of zero ([3]). 

3 Results and Discussion 

In the previous section, the general idea as well as the detailed steps for the 
construction of the software sensor have been given. To complete the anal- 
ysis of the approach, this section presents the experimental results obtained 
at the laboratory level at which the software sensor was tested. 

In total, six series of experiments of Mezcal fermentation at temperature 
of 32°C were performed. The temperature has been maintained constant 
using a heat exchanger device, see Fig. [TJ First, we performed three experi- 
ments to generate the necessary biomass and ethanol concentration data for 
the training of the ANN. Once the error given by the ANN goes below the 
value of 0.003 biomass units, the adaptive scheme was added to complete 
the software sensor procedure. Then, three more experiments were carried 
out for further testing the scheme. The results concerned with the perfor- 
mance of the ANN to infer the ethanol and biomass concentrations from 
the redox potential measurements are presented in Figs. [2] and [3j where one 
can appreciate the degree of accuracy of the algorithm. Notice that Fig. 
[3] presents several slope changes that do not allow a direct inference of the 
presence of a mixture of microorganisms during the fermentation process. 

Although the main goal is to estimate n(t), the biomass estimation is 
also performed, and the difference between this estimation and the value pre- 
dicted by the ANN procedure is used as a correction term in the algorithm. 




(3) 



S 




time (hr) time (hr) 



Figure 2: The right plot illustrates the time evolution of the redox potential 
in the Mezcal fermentation process. In the left plot the blue solid line rep- 
resents the ANN-estimated ethanol concentration while the red dots stand 
for the experimental values. 

Figure [5] shows the ability of the software sensor to rebuild the biomass con- 
centration data. It should be highlighted that the simplicity associated with 
the implementation of this algorithm and the fact that a single measured 
signal, the redox potential, is needed are very promising features from the 
technological and industrial standpoints. 

In Fig. [5j it is possible to distinguish three regions, labeled as A, B and 
C. In general terms, in region A one is not able to get any conclusion on 
the dynamical behavior of the process since during the corresponding time 
span the observer has a transient behavior. The span of time corresponding 
to region A is given by the time that the observer takes up to minimize the 
error between its estimated biomass concentration and that given by the 
ANN procedure. 

On the other hand, we associate regions B and C with the presence of 
two different groups of microorganisms, namely bacteria and yeasts. This 
assumption is based on the fact that the growth rate of bacteria is faster than 
that of yeasts. Thus, from Fig. the presence of a mixed culture growth 
comes out naturally, a fact which is not so obvious when one examines only 
the biomass data given in Fig. [3j Note in addition that the end of the 
fermentation process is quite clear in Fig. 
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Figure 3: The continuous lines correspond to the estimated biomass data 
and the dots to the experimental data. Blue, red and green colors correspond 
to three of the experiments carried out in our laboratory. 

4 Conclusions 

In this work, we have shown that in the Mezcal fermentation process the 
redox potential can give relevant information on the metabolic activity of the 
involved microorganisms, including both ethanol and biomass concentration. 
Besides, the problem of estimating the specific growth rates in the Mezcal 
fermentation process has been addressed through a combination of a simple 
ANN procedure and a nonlinear observer algorithm. Moreover, the strategy 
presented here detects in a clear way the end of this fermentative process. 
The latter fact is quite relevant from the production point of view since 
it is the basic piece of information by which one can save time and avoid 
the degrading process due to the conversion of the ethanol into acetic acid. 
This leads to an improvement in the quality of the product. In addition, 
the computational scheme provides a very appropriate tool for a fast and 
reliable quality control and can be used to ensure the homogeneity of the 
final product. The methodology presented in this paper is general and can 
be also used in automatic control applications. 



10 



q| , , 1 q| , , 1 q| , , 1 

5 10 5 10 5 10 

time (hr) 

Figure 4: (Blue) solid lines represent the estimated biomass concentration 
given by the software sensor and (red) dash-dotted lines stand for the ANN- 
predicted biomass concentrations. From left to right, the displayed plots 
correspond to the three experiment performed to test our approach. 
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